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ORIGINAL ARTICLE 

Transcripts involved in calcium signaling and telencephalic 
neuronal fate are altered in induced pluripotent stem cells from 
bipolar disorder patients 

HM Chen 1 ' 5 , CJ DeLong 2 ' 5 , M Bame 1 , I Rajapakse 3 , TJ Herron 4 , MG Mclnnis 1 ' 6 and KS O'Shea 2 ' 6 

Bipolar disorder (BP) is a chronic psychiatric condition characterized by dynamic, pathological mood fluctuations from mania to 
depression. To date, a major challenge in studying human neuropsychiatric conditions such as BP has been limited access to viable 
central nervous system tissue to examine disease progression. Patient-derived induced pluripotent stem cells (iPSCs) now offer an 
opportunity to analyze the full compliment of neural tissues and the prospect of identifying novel disease mechanisms. We have 
examined changes in gene expression as iPSC derived from well-characterized patients differentiate into neurons; there was little 
difference in the transcriptome of iPSC, but BP neurons were significantly different than controls in their transcriptional profile. 
Expression of transcripts for membrane bound receptors and ion channels was significantly increased in BP-derived neurons 
compared with controls, and we found that lithium pretreatment of BP neurons significantly altered their calcium transient and 
wave amplitude. The expression of transcription factors involved in the specification of telencephalic neuronal identity was also 
altered. Control neurons expressed transcripts that confer dorsal telencephalic fate, whereas BP neurons expressed genes involved 
in the differentiation of ventral (medial ganglionic eminence) regions. Cells were responsive to dorsal/ventral patterning cues, as 
addition of the Hedgehog (ventral) pathway activator purmorphamine or a dorsalizing agent (lithium) stimulated expression of 
NKX2-1 (ventral identity) or EMX2 (dorsal) in both groups. Cell-based models should have a significant impact on our understanding 
of the genesis and therefore treatment of BP; the iPSC cell lines themselves provide an important resource for comparison with 
other neurodevelopmental disorders. 
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INTRODUCTION 

Bipolar disorder (BP) is a severe neuropsychiatric disease 
characterized by recurrent episodes of manias, pathologically 
energized states with misguided volition and behavior in a mood 
state of intoxicating euphoria (or irritability) and depressions that 
are low moods with compromised energy, volitional states and 
diminished cognitive capacity. The social, personal and vocational 
consequences are frequently disastrous, 1 and the WHO ranks BP 
among the top causes of lifelong disability. 2 

The heritability of BP is estimated at 80% and adoption, twin 
and family studies have consistently identified a significant 
genetic contribution; 3,4 many genes 5 and genetic mechanisms 6 
have been implicated in BP, consistent with multiple disease loci, 
each contributing modest effects. 7,8 Genome-wide association 
studies have identified a few replicated loci including polymorph- 
isms in CACNA1C, the alpha subunit of the L-type voltage-gated 
calcium channel, the strongest and most replicated risk. 5,9 " 13 
Additional replicated polymorphisms have been identified in 
ANK3, which regulates assembly of voltage-gated sodium 
channels at the node of Ranvier, 14,15 and SYNE1, which has a 



critical role in synaptic clustering. 16 Overall, available data suggest 
that alterations in membrane organization and/or ion channelo- 
pathies may contribute to BP; in fact, dysregulation of Ca 2+ 
signaling is one of the most reproducible observations in BP. 17 

Like many complex central nervous system (CNS) disorders, the 
roles of neuropathology, environment and trauma in BP are 
largely unknown; there are reports of ventricular enlargement, 
gray matter, hippocampal and cerebellar hypoplasia, and anoma- 
lies of white matter tracts, 18-21 that are nonspecific and 
inconsistent. There is also evidence of abnormal neuronal cell 
migration 22 and of altered adult neurogenesis. 23,24 Because viable 
brain tissue is unavailable, and there is no well-defined genetic 
event or neuropathology, there are simply no good neural models 
to study BP. 

Since subtle alterations in developmental pathways can 
produce neurological consequences that only become apparent 
much later in life, 25,26 a developmental model to study neuronal 
differentiation in BP is required. Because they have tri-lineage 
differentiation potential similar to human embryonic stem cells 
(hESCs), induced pluripotent stem cells (iPSCs) provide the 
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opportunity to study the stepwise differentiation of patient- 
derived cells into neurons and glia to identify alterations in cell 
behavior and test novel therapeutic approaches. 27,28 In BP, they 
offer the unique opportunity to obtain a robust source of cells 
from individuals with prospective longitudinal observations, 
individuals on or off medications, patients who have had multiple 
episodes, different diagnoses and treatments, unlike hESC. iPSCs 
have been derived from patients with neurodegenerative 
disorders 29-34 and other neuropsychiatric disorders, 35-38 but are 
just being developed for BP. 

The current investigation demonstrates that iPSC from BP 
patients are similar to controls (C) in their tri-lineage differentia- 
tion capacity and gene expression profiles; however, the expres- 
sion of membrane receptors and ion channel transcripts was 
significantly increased in BP neurons, and different sets of 
neuronal lineage-specific transcription factors were induced. 
Consistent with changes in membrane composition, exposure of 
BP neurons to lithium produced significant alterations in wave 
amplitude and calcium transient compared with control neurons. 
Transcripts involved in dorsal telencephalic fate (the default 
identity) were enriched in control neurons, whereas there was an 
increase in genes involved in ventral fate and GABAergic 
interneuron differentiation in BP neurons. The dorsal/ventral 
neuronal identity of both groups could be influenced by exposure 
of iPSC to dorsalizing or ventralizing factors. This work provides 
the first information on gene expression during neuronal 
differentiation in BP, a rich resource to identify alterations in the 
differentiation of glia and neurons, as well as an opportunity for 
comparisons with other neuropsychiatric disorders. 



MATERIALS AND METHODS 

Patient population 

The clinical sample was obtained from the Prechter Longitudinal Study of 
Bipolar Disorder. 39 Subjects contributing a skin sample were from a 
psychiatric clinic in a mid-western college city, Caucasian, and were 
diagnosed with Bipolar 1 Disorder, or healthy unaffected controls were 
ascertained through advertising on the University of Michigan Clinical 
Studies website. All participants were assessed clinically using the 
Diagnostic Interview for Genetic Studies. 40 BP patients were determined 
by clinical evaluation to have a chronic disorder with mood symptoms 
most of the time, age of onset ranging from 15 to 23 years. Two of three 
Bipolar 1 Disorder patients were considered lithium responsive and all had 
exposures to antidepressant, anticonvulsants and antipsychotic medica- 
tions. One patient had comorbid panic disorder and a second patient was 
treated for metabolic syndrome with a statin (rosuvastatin) and an oral 
hypoglycemic (pioglitazone). A second patient was receiving levothyrox- 
ine. Two patients and two controls were 50 ± 5 years age and the other pair 
were aged 30 and 34 years at time of sample. Four patients were male and 
two were female. 

Fibroblast derivation 

Skin biopsies were obtained in the UM Clinical Research Unit. Samples 
were dissociated in trypsin, expanded in fibroblast medium comprising 
Dulbecco's Modified Eagle Medium (DMEM, 11965-092; Invitrogen, Grand 
Island, NY, USA) containing 10% heat-inactivated fetal bovine serum 
(Invitrogen, 10082-147), 1% non-essential amino acids (NEAA) (Invitrogen, 
11140-050), 2mivi glutamax (Invitrogen, 35050-061) and 1x Penicillin/ 
Streptomycin solution (Invitrogen, 15140-122) to obtain dermal fibro- 
blasts. Stock fibroblasts from each patient were banked and cells were 
analyzed for their expression of fibroblast-restricted genes, pluripotency 
factors, lineage markers and mycoplasma. 

Virus preparation 

Retrovirus was made by transfecting 293FT cells with pUMVC and pCMV- 
VSV-g plasmids, and individual plasmids: pMXs-hKLF4, pMXs-hSOX2, pMXs- 
hOC74 or pMXs-hcMYC (Addgene, Cambridge, MA, USA) at a 1:1:1 ratio 
using Lipofectamine 2000 (Invitrogen, 11668-019). Virus supernatant was 
collected and combined, filtered through a 0.45-um filter (EMD Millipore, 
SCHVU01 RE, Billerica, MA, USA) and concentrated 80-100-fold in a 



centrifugal filter with a 1 00 000 MWCO (EMD Millipore, UFC910008). The 
concentrate containing all four factors was then used immediately to 
transduce fibroblasts at a final concentration between 3- and 12-fold. 

iPSC derivation 

Fibroblasts at passage 5-6 were plated at 5x10 4 cells in 35-mm tissue 
culture dishes in fibroblast medium and transduced with retroviral 
constructs expressing the pluripotency factors OCT4, SOX2, KLF4 and 
c-MVC 41 Fibroblasts were transduced on day 1 after plating, with daily 
medium changes from day 3 to day 6. On day 7, forming iPSC were 
trypsinized and passaged to mouse embryonic fibroblast (MEF) feeders for 
expansion in hESC medium compromising DMEM/F12 (Invitrogen, 11320- 
-032), 20% KOSR (Invitrogen, 10828-028), 4 ng ml" 1 FGF2 (EMD Millipore 
GF003), 1% NEAA, 2 mivi glutamax and 0.1 mivi p-mercaptoethanol (Sigma- 
Aldrich, St Louis, MO, USA, M3148). After 10 days on MEFs, medium was 
changed to MEF-conditioned hESC medium, then hESC-like colonies 
manually dissected and passaged to new MEFs for expansion in hESC 
medium. From each fibroblast line, 5-10 colonies were selected for 
expansion to iPSC lines; a minimum of 5-6 iPSC lines were derived from 
each fibroblast (patient) sample. 

Pluripotency characterization 

Analysis of the expression of pluripotency factors NANOG, 50X2 and OCT4, 
and cell surface expression of SSEA3, SSEA4, TRA-1-60 and TRA-1-81, and 
alkaline phosphatase was carried out in immunocytochemistry and 
polymerase chain reaction (PCR) (Supplementary Figure 1). Controls for 
pluripotency and differentiation capability include parallel analysis of a 
well-characterized iPSC line (UM-C2) and H7 or UM4-6 hESC lines. 

Differentiation capacity 

Assessment of the tri-lineage differentiation capability of the iPSC lines was 
done by 45 days of culture in non-adherent dishes in hESC medium lacking 
FGF2 to form embryoid bodies (EBs). Lineage differentiation was assessed 
using cell-type restricted antibodies and PCR (Supplementary Figure ID, 
Supplementary Tables 1 A and B). Ectoderm was identified using nestin or 
(3-III tubulin; mesoderm: Brachyury, smooth muscle actin or decorin; and 
endoderm: Sox17, alpha-fetoprotein or Foxa2. 

Neuronal differentiation 

iPSCs were manually dissected from MEF, then grown in suspension 
culture for 4 days in EB medium supplemented with 2 um dorsomorphin 
(R&D Systems, 3096; Minneapolis, MN, USA) and 5 mm SB431542 (R&D 
Systems, 1614). Medium was then changed to neural induction medium 
consisting of DMEM/F12, 1 x N2 supplement (Invitrogen, 17502-048), 1 x 
NEAA, 2 itim glutamax, 0.1 mivi (3-mercaptoethanol, 2 um dorsomorphin, 5 
um SB431542 and 20 ng ml -1 FGF2 for an additional 3 days in suspension. 
Neurospheres were then planted on laminin (Sigma-Aldrich, L2020)-coated 
culture dishes for rosette formation. Rosettes were manually picked 2-3 
times, dissociated in Accutase then re-plated on poly-D-ornithine-laminin- 
coated substrates for neuronal differentiation in medium composed of 
neural basal medium (Invitrogen, 21103-049), 1 x B27 supplement 
(Invitrogen, 17504-044), Iitim glutamine, 1% NEAA and 1x Penicillin/ 
Streptomycin for an additional 4, 8 or 12 weeks. Medium was changed 
every other day. To determine whether dorsal-ventral fate could be 
respecified, cells were grown with the Smoothened agonist (Hedgehog 
pathway activator) purmorphamine (1 um; Cayman Chemical, Ann Arbor, 
Ml, USA, 10009634), or with the dorsalizing agent lithium chloride (LiCI; 
1 itim; Sigma-Aldrich). Controls were exposed to DMSO (carrier) alone. 42,43 

Reverse transcriptase-polymerase chain reaction 
Total RNA from iPSCs, EBs or neurons was extracted with Trizol (Invitrogen, 
15596-0926), ±5ug of UltraPure glycogen (Thermo Fisher Scientific, 
Waltham, MA, USA) to enhance RNA recovery. RNA was treated with 
DNase (Thermo Fisher Scientific) and reverse transcribed (Superscript III, 
Invitrogen) using random nonamers. The cDNA was then used in PCR 
(primers in Supplementary Table IB). 

Microarray analysis 

To characterize the iPSC and to determine whether there were differences 
in their gene expression profiles with neuronal differentiation, total RNA 
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Figure 1. Derivation, characterization and neuronal differentiation of iPSC. Fibroblasts (a) express Te-7 (b; Cy3-secondary antibody, red). 
Rep rog rammed colonies of iPSC form and are passaged to MEFs. They downregulate Te-7 (c; Cy3-secondary antibody), and exhibit nuclear 
expression of Nanog (d; FITC-secondary antibody). For differentiation, cells are removed from MEF and grown in defined medium as 
neurospheres (e), then plated to form rosettes (f). Rosettes are picked twice, then differentiated in adherent culture in defined medium 
forming nestin + (g; FITC-secondary antibody) neural progenitor cells, and with additional time in vitro undergo widespread differentiation to 
p-lll tubulin + neurons (h; Cy3-secondary antibody), with punctate expression of Synapsin 1 (Syn1) along their processes (i, Cy3-secondary 
antibody). All cells, but a, are stained with Hoechst 23487 (blue) to identify nuclei. Scale bars in a-d,g = 500 mm; e = 2 mm, f,h,i = 200 mm. 



was isolated from six individual iPSC cell lines (3 BP patients and 3 controls 
(Q) before and following 8 weeks of neuronal differentiation using the 
TRIzol reagent (Invitrogen, 15596-026). RNAs were DNased followed by the 
RNeasy MinElute Cleanup reagent (Qiagen, Hilden, Germany, 74204). RNA 
quantity and quantity were measured in an Agilent 2100 bioanalyzer 
(Agilent Technologies, Santa Clara, CA, USA) producing RNAs with an 
integrity measure > 9.5. They were amplified and hybridized to GeneChip 
U133 Plus 2.0 microarrays (Affymetrix, Santa Clara, CA, USA). Data quality 
control parameters included scaling factor, noise, background, percentage 
of present calls, 573' signal ratios observed in GAPDH mRNA and 573' 
signal ratios of spiking genes. Raw data (U133P2 CEL files) were processed 
using the Bioconductor package 'Affy' for R. Normalized Log 2 transformed 
probe intensity Robust Multi-chip Average values were used in differential 
gene expression analyses and transcripts with significantly altered 
expression (P^0.05) between BP and controls were identified using 
f-test with R and David software. 44 To assess similarities between groups, 
the Symmetrized Stein Distance was calculated between estimated sample 
covariance matrices for each sample pair (BP vs control iPSC, BP vs control 
neurons and iPSC vs neurons). Statistical significance was calculated using 



the permutation test. 45 Validation of the microarray results was performed 
with TaqMan 'Assays on Demand' (ABI, Foster City, CA, USA) for seven 
genes selected based on fold changes between BP and control samples, 
using the 2~ AACt method. 46 Only complete sets of iPSC (3 BP and 3 C 
patients) and neurons (same 3 BP and same 3 C cell lines from 6 
individuals) were analyzed to minimize stochastic changes due to culture 
conditions. 



Calcium signaling 

To simultaneously measure action potential and wave propagation, we 
differentiated neurons as described above for an additional 4-8 weeks, 
loaded them with the intracellular calcium-sensitive dye Fluo-4AM, (10 mivi, 
Molecular Probes, Eugene, OR, USA). After 30 min, non-bound dye was 
removed and cells were examined in a Nikon A1R laser scanning confocal 
microscope system (Nikon Corporation, Chiyoda, Tokyo, Japan) to identify 
spontaneous activity. To evoke action potentials, 50 mivi KCI was added to 
the dish and signaling was recorded. 47 To determine whether lithium 
pretreatment would alter signaling, sets of neurons (3 BP and 3 C) were 
cultured in the presence of LiCI (1 itim) for 24 h before loading and image 
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acquisition. For image analysis, regions of interest were defined as cell 
bodies or axons, and wave amplitude and propagation were measured 
from 5 to 10 regions per sample. 



RESULTS 

Cells 

Fibroblast lines have been derived from 22 BP and 10 C patients, 
expanded and banked to 12 stock vials (or cryostraws) from each 
patient. From 3 BP and 3 C fibroblast cell lines, we have 
reprogrammed at least 3 iPSC lines per individual (BP iPSC = 28; 
C iPSC=14). With reprogramming, fibroblast-restricted genes 
were downregulated (Figures 1a-c; Supplementary Figure 1A) 
and expression of pluripotency factors was induced (Figure 1d; 
Supplementary Figures 1B and C). Cells expressing markers of all 
three germ layers (ectoderm, endoderm and mesoderm) formed 
following EB differentiation (Supplementary Figure 1D). With 
differentiation of neurospheres (Figure 1e) followed by neural 
rosettes (Figure 1f) Sox1, Sox2, Sox3, Pax6 and nestin were 
expressed (Figure 1g, Supplementary Figure 2). After 4 weeks, 
neurons with dendrites and axons expressing MAP2, TauAB and 
P-lll tubulin (Figure 1h) formed. Neurons adopted a largely 
forebrain phenotype, with VGIut+ pyramidal-shaped neurons, 
GABA+ interneurons, scattered serotonergic and dopaminergic 
neurons present in the cultures. Many expressed punctate 
Synapsinl along their neurites (Figure 1i), with little glial fibrillary 
acidic protein (GFAP) in the cultures. 

Microarray analysis 

Pooled data: BP and control. Analysis of the 54 676 probesets 
identified 20 550 unique transcripts expressed in iPSC and 20 718 
in neuronally differentiated cells. We replicated five out of seven 
genes selected from our microarray in quantitative reverse 
transcriptase-polymerase chain reaction, supporting the validity 
of these results (Supplementary Table 1C). In the iPSC samples, 
2538 transcripts were expressed at significantly different levels 
(Supplementary Figure 3A). KEGG pathway analysis identified one 
significant pathway in the iPSC, calcium signaling, n = 39, 
P^2.5x10~ 2 , Bonferroni correction. Many of these transcripts 
are associated with neuroplasticity and development, others with 
stress-activated pathways (Supplementary Table 2). Of the 20 718 
transcripts identified in neuronally differentiated iPSC, 1835 
transcripts were differentially expressed between groups (nominal 
P^0.05), but KEGG pathway analysis did not identify significantly 
altered pathways in the neuronal population. Additional analysis 
of the differentially expressed transcripts did not identify 
significant differences between BP vs control iPSC (distance = 
0.036, 0.072), not surprisingly pooled iPSC vs neurons were 
significantly different (distance = 0.1 74, P^ 0.01 2). Interestingly, BP 



and control neurons were also significantly different in their 
transcriptomes (distance = 0.1 03, 0.041); Supplementary Figure 
3B. 

Within-groups analysis: BP iPSC vs BP neurons and control iPSC vs 
control neurons 

Transcripts expressed in iPSC and neurons: There were 14 993 
unique transcripts that were significantly differentially expressed 
in BP iPSC compared with neurons derived from them; there were 
15 311 transcripts in control iPSC vs control neurons that were 
significantly altered. Of these, 7242 transcripts were expressed at 
twofold or higher levels in BP iPSC, 5120 were expressed at ^ 
twofold levels in BP-derived neurons. In control cells, 7098 unique 
transcripts were expressed at twofold or higher levels in iPSC and 
5089 at ^ twofold levels in neurons. 

i. Transcripts enriched in iPSC Functional Annotation Chart 
Analysis identified 475 Gene Ontology (GO) terms associated with 
genes expressed at significantly higher levels in BP iPSC than in 
differentiated BP neurons; 463 terms in control iPSC vs neurons. 
These included transcripts involved in cell cycle control, RNA 
binding and genes involved in development (Table 1 summarizes 
clusters involved in early differentiation). 

Both sets of iPSC express similar pluripotency factors, positive 
regulators of cell cycle, genes involved in very early embryonic 
development, signaling molecules, and factors involved in 
proliferation. They express few anterior-posterior patterning 
factors (HOX) and numerous genes involved in chromatin 
modification. There were slight variations in the expression of 
individual transcripts between groups, but the functional sig- 
nificance, if any, of these differences is not known. 

ii. Transcripts enriched in neurons Transcripts associated with 425 
Gene Ontology categories were significantly increased in BP 
neurons compared with BP iPSC; 437 in control neurons compared 
with control iPSC (Table 2 describes categories involved in 
neuronal differentiation). 

With neuronal differentiation, transcripts associated with 
pluripotency, positive control of the cell cycle, and many 
chromatin remodeling factors were downregulated, whereas 
genes involved in signaling, transcripts associated with anterior- 
posterior patterning (HOX genes) and SOX genes were induced 
with differentiation. Cell surface molecules, many involved in 
neurite outgrowth, and transcripts that encode membrane 
channels were also upregulated with neuronal differentiation. 
Transcripts involved in early neuronal differentiation including 
ELAVL1, 2,3,4) FOXP1; MAP2; MAPT; MBD2; NEFL, NEFM and NRG1-4 
were increased, whereas genes associated with glial lineage 
differentiation were not identified in BP, but GFAP was present in 
controls. Transcripts previously suggested to be dysregulated in 
BP, including ANK3, CACNA1C, GSK3/5 (glycogen synthase kinase 3 



Table 2. Enriched in neurons 






BP neurons 




Control neurons 


GO term 


n 


P*^ 




n P**c 


Synapse 


183 


2.1 xlO 


-22 


180 1.5x10~ 21 


Neuron projection 


173 


3.4x10 


-20 


176 1.2X10" 21 


Neuron development 


161 


4.5x10 


-12 


167 7.2X10" 15 


Voltage-gated cation 


66 


7.8x10 


-4 


61 1.6x10" 2 


channel activity 











Abbreviation: BP, bipolar disorder; GO, Gene Ontology; * Bonferroni 
correction. Neurons differentiated from BP and control iPSC expressed 
similar sets of transcripts, many associated with neuronal differentiation, 
synapses and channel activity. 



Table 1. Enriched in iPSC 



BP iPSC Control iPSC 



GO term n P* ^ n P* ^ 



Blastocyst development 32 5.1 x10~ 4 27 2.1 x10 -3 

Post-transcriptional regulation 109 5.1 x10~ 3 109 2.3x10~ 3 
of gene expression 

In utero development 90 2.5x10~ 2 85 4.2x10~ 3 

Abbreviations: BP, bipolar disorder; GO, Gene Ontology; iPSC, induced 
pluripotent stem cells; *, Bonferroni correction. iPSC from BP and control 
individuals expressed very similar sets of transcripts compared with 
neurons differentiated from them. Functional Chart analysis identified 
comparable sets of GO terms and similar transcripts in each category, 
particularly genes associated with early embryonic development. 
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Table 4. Analysis of neurons 


//'/'. Enriched in BP neurons 


/V. Enriched in C 








neurons 


GO term 


n 


P*^ 


Go term n P*^ 


Plasma 


38 


4.5x10~ 3 


Nucleus 44 5.6x10~ 2 


membrane part 








Ion transport 


15 


1.3x10 -2 





Abbreviations: BP, bipolar disorder; C, control; GO, Gene Ontology; 
*, Bonferroni correction, iii: Functional Annotation Cluster analysis 
identified two significant clusters of transcripts expressed at higher levels 
in BP vs control neurons, iv: There was a single cluster containing nuclear 
transcripts that were expressed at higher levels in control neurons. 



beta), MECP2, NCAN, SYNE1 and ZNF804A, were significantly 
upregulated as both control and BP iPSC differentiated to neurons, 
whereas BDNF and DISCI were not changed. 

The ribonuclease DICER1, which is involved in microRNA 
processing, was significantly higher in BP but not control neurons. 
Factors involved in ventral neuronal fate such as NKX2-1 and PAX2 
were significantly induced in BP, whereas transcripts involved in 
dorsal fate such as EMX2 (dorsal telencephalon), and HOX genes 
{HOXA 1,6,7) required for positional identity and regional differ- 
entiation of the hindbrain were significantly higher in control 
neurons. 

These data demonstrate that with neuronal differentiation, cells 
withdraw from cell cycle, express different sets of signaling 
molecules than iPSC and express largely overlapping sets of 
neuronal lineage markers. Enrichment of dorsal telencephalic 
transcripts, such as PAX6 and EMX2, and the cortical determinate 
LHX2 (ref.48) in control and expression of ventral determinates, 
such as PAX2 and NKX2-1 (basal forebrain progenitors 49 ), involved 
in formation of the medial ganglionic eminence and GABAergic 
interneuron differentiation 50 in BP neurons suggest that neuronal 
identity is altered, as mutual inhibition between PAX6 and NKX2-1 
determines dorsal/ventral boundaries in the telencephalon. 51 

Between-groups comparisons: BP vs control cells 
Analysis of pluripotent cells: There were 2547 unique transcripts 
expressed at significantly different levels in iPSC derived from BP 
vs control individuals. 

i. Transcripts enriched in BP vs control iPSC One-hundred eighty- 
six transcripts were expressed at 1.5-fold or greater levels in iPSCs 
derived from individuals with BP compared with control iPSCs. 

There were 14 significant clusters of genes (Table 3i) in this 
group. Many transcripts previously associated with neurodevelop- 
mental/neuropsychiatric disorders or susceptibility loci for BP were 
identified (Supplementary Tables 3-7). 

ii. Transcripts enriched in control vs BP iPSC Forty-six transcripts 
were expressed 1.5-fold or higher in iPSC from control vs BP 
individuals, including genes involved in cell polarity and in lineage 
differentiation, including the TOR complex2 factor RICTOR, which 
regulates cell growth. There was one cluster of transforming 
growth factor-(3-related transcripts that was significantly increased 
(Table 3ii). Members of this cluster included GDF3 and INHBE that 



Table 3. Analysis of iPSCs 



have been associated with proliferation and pluripotency and are 
expressed in the ICM and/or in hESC, 52 whereas LEFTY1 and 
LEFTY2 are involved in nodal inhibition, neural induction and early 
embryo patterning. 53 None was previously associated with 
psychiatric or human neurodevelopmental disorders. These data 
indicate that iPSC from BP and control cells are strikingly similar in 
their transcriptional repertoire. No genes associated with the core 
regulatory mechanisms that maintain pluripotency were differen- 
tially expressed in these two populations. 

iii. Genes enriched in BP vs control neurons Of the 1312 unique 
transcripts that were altered in neurons, 140 transcripts were 
increased 1.5-fold or greater in BP, including transcripts associated 
with cell-cell and/or cell-matrix interactions and with axon 
outgrowth, with neuronal differentiation, organization of 
synapses, secreted growth factors, and neurotransmitters and 
their receptors. Many membrane-associated channel genes, 
transcripts involved in methylation and control of gene expression 
and with the cytoskeleton were also increased in BP neurons. 
Functional Annotation Cluster Analysis identified plasma mem- 
brane part and ion transport as significantly increased in BP 
neurons (Table 4iii). Both clusters (Supplementary Tables 8 
and 9) contained transcripts previously associated with neurode- 
velopmental/neuropsychiatric disorders, as well as BP 
susceptibility loci. 



/'. Increased in BP iPSC ii. Increased in C iPSC 



GO term 


n 


P**C 


GO term n P* ^ 


Embryonic morphogenesis 


17 


9.0X10" 5 


TGFp related 4 5.7x10~ 4 


Signal 


57 


6.1 x10~ 4 




Axonogenesis 


10 


9.4X10" 3 




Appendage development 


8 


6.7X10" 3 




Skeletal system development 


12 


1.7X10" 2 




Cell adhesion 


20 


6.7X10" 3 




Female gonad development 


7 


6.1 x10~ 3 




Skeletal system morphogenesis 


7 


3.0X10" 2 




Positive regulation of developmental process 


13 


4.9 xlO" 3 




Cell-cell adhesion 


10 


4.9x10~ 2 




Transcription from Polll 


18 


3.0X10" 2 




Blood vessel development 


10 


2.8X10" 2 




Transcription regulation 


36 


2.2X10" 2 




Tube development 


9 


4.5X10" 2 




Abbreviations: BP, bipolar disorder; C, control; GO, Gene Ontology; iPSC, induced pluripotent stem cells; TGFp, transforming growth factor-p; * Bonferroni 
correction, i: Functional Annotation Cluster analysis identified 14 clusters of transcripts that were expressed at significantly higher levels in BP vs control iPSC, 
most associated with early development, ii: A single cluster of transforming growth factor-p-related transcripts was expressed at significantly higher levels in 



control iPSC. 
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iv. Enriched in control vs BP neurons Ninety-five transcripts were 
expressed 1.5-fold or higher in control vs BP neurons, including 
genes involved in cytoskeletal organization and function, cell 
surface/extracellular matrix factors, apoptosis, DNA repair, in cell 
cycle, and in signaling. Genes involved in early neuronal 
differentiation and in dorsal patterning of the CNS including 
EMX2, EOMES, TCF3 as well as VGLUT1 (SLC17A7; previously 
reported to be decreased in BP 54 ) were enriched in control but 
not BP neurons; there was a single significant cluster of nuclear 
factors (Table 4iv; Supplementary Table 10). 

Several patterning factors were altered in BP vs control neurons. 
HOXA1, A6 and A7 that are expressed in the cerebellum and 
brainstem were significantly increased in control neurons, 
consistent with clinical reports of cerebellar hypoplasia 20,21,55 in 
BP. Of transcripts associated with circadian rhythm, TIPIN, PENK 
and TAC1 were expressed at significantly higher levels in BP vs 
control neurons; whereas in iPSC NPAS2 was expressed at 
significantly lower levels in BP vs controls; CRY! and RORA were 
expressed at higher levels. 



Comparisons with other data sets 

We previously examined gene expression in postmortem RNA 
samples from the premotor cortex of controls, from individuals 
with BP taking typical antipsychotic medications (BP+), and from 
unmedicated patients with BP (BP-). Remarkably, gene expression 
profiles of BP+ were more similar to controls than to BP-. 56 When 
differentially expressed transcripts from BP- brain were compared 
with those from BP iPSC-derived neurons, 151 significantly altered 
genes were identified; 57 genes were expressed in a similar 
direction (x 2 = 651; P^1x10~ 7 ) including transcripts involved in 
cell migration, calcium signaling and neurotransmitter receptors 
(Supplementary Table 11). 

Compared with the genes enriched in schizophrenia patient- 
derived iPSC neurons, 35 nine were significantly different in bipolar 
vs control neurons in our study. 50X6 and FGFR3 were down- 
regulated in both, whereas FGF14, LGALS8, PDE4D, RHOBTB3, 
SLC41A2, UNC5C and ZMAT4 were expressed at higher levels in 
schizophrenia and BP than in control neurons. Compared with an 
analysis of neurons differentiated from iPSC from patients with 
Timothy Syndrome {CACNA1C mutation 37 ), 233 genes were 
expressed in a similar pattern as in BP neurons, including many 
transcripts involved in neuronal differentiation. Clearly, it will be 
important to compare our analyses with those of additional BP 
tissue, cell and mouse models. 



Patterning 

The MGE marker NKX2-1 was expressed within rosettes and in 
cells rapidly migrating from them, whereas EMX2 was expressed 
at the base of the rosette where cells were beginning to 
differentiate (Supplementary Figure 4). Control neurons expressed 
more EMX2 and less NKX2-1 than BP neurons regardless of culture 
conditions; PAX6 was expressed in rosettes in both groups. Both 
BP and control iPSC could be patterned by exposure to the 
Hedgehog pathway activator purmorphamine by increasing the 
number of NKX2-1 -positive cells (Supplementary Figures 4A and 
D) or to LiCI by increasing the number of neurons expressing 
EMX2 (Supplementary Figures 4C and F). 



Calcium signaling 

By 12 weeks, spontaneous neuronal firing could be observed in 
these cultures. Lithium pretreatment of BP neurons significantly 
decreased their calcium transient (101.65 ±31.64 vs 290.22 
±0.52 sec; 0.007) and wave amplitude (0.76 ±0.07 vs 
1.89 ±0.26 AF/F 0 ; 0.004) compared with lithium-exposed 
control neurons (Supplementary Figure 5). We also observed that 
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clusters of soma were capable of clearing Ca2 + more rapidly than 
axon- and dendrite-rich regions. 



DISCUSSION 

There is a clear need for cell-based models of BP to test genetic 
and epigenetic factors, the role of the environment (particularly 
stress), to determine when neuronal differentiation is affected, and 
to identify and test new therapeutic approaches. The recent 
discovery that four transcription factors can reprogram adult 
somatic cells into PSCs 41 has provided an important tool to study 
human disease. 31,57,58 In addition to interrogating the many 
hypotheses regarding the etiology of BP, the iPSC will be useful in 
modeling and comparing with other neuropsychiatric disorders 
that have developmental or affective phenotypes. The banked 
fibroblast lines also provide a valuable resource for direct 
reprogramming of the desired cell or neuronal subtype and for 
alternative reprogramming technologies. 

The expression of factors associated with pluripotency as well as 
their tri-lineage differentiation capability, support the successful 
reprogramming of dermal fibroblasts to iPSC. Not surprisingly, BP 
and control iPSCs did not differ significantly in their transcriptional 
profiles. Although there were many similarities in gene expression, 
there were differences that will allow us to begin to examine the 
pathways previously suggested to be affected in BP and those 
identified in this study including calcium signaling and microRNA 
processing. 

Neurons differentiated from BP iPSC were significantly different 
in their gene expression profile from those derived from control 
iPSC. BP neurons expressed more membrane receptors and ion 
channel genes than control neurons, particularly transcripts 
involved in calcium signaling. Subtle dysregulation of calcium 
signaling has been suggested to generate inappropriate neuronal 
responses and enhance the tonic excitability that maintains brain 
rhythms in BP. 59 Calcium signaling also has a critical role in 
development, synaptic plasticity and homeostasis of the CNS; 
misregulation during development can produce subtle but 
widespread alterations in lineage differentiation and plasticity 
throughout the nervous system and may influence susceptibility 
to BP. 60 Our preliminary analyses suggest that calcium signaling is 
also altered in BP neurons and is sensitive to lithium pretreatment. 
Lithium acts both pre- and post-synaptically, can alter calcium 
uptake and intracellular levels, and via G5/G/3 can be neuropro- 
tective. It will now be possible to examine the detailed effects of 
lithium exposure as well as effects of other BP drugs and other 
signaling pathways in these cells. 

These results support the concept that many complex, 
neuropsychiatric disorders that are identified later in life, including 
gp 61-63 ma y resu | t f rom alterations in neural differentiation that 
occur during development. BP neurons expressed transcripts 
involved in the differentiation of ventral neuronal subtypes, 
including FOXP2 and NKX2-1, the hallmark of the medial 
ganglionic eminence, 51 whereas controls expressed higher levels 
of factors that confer (or maintain) dorsal telencephalic neuronal 
identity, including EMX2, PAX6, TBR2, TCF3 and ZNF536. As 
telencephalic neurons are involved in a number of neurodegen- 
erative (Huntington's disease and Alzheimer's disease) and 
neurodevelopmental conditions (autism spectrum disorder, schi- 
zophrenia, BP and Rett syndrome), alterations in subtype 
specification may begin to explain the origin or susceptibility to 
these conditions. 64 Because lithium activates WNT pathway 
signaling, dorsalizing early CNS progenitors, 26 lithium treatment 
may re-direct neural stem cell fate to dorsal cortical derivatives. 
Many additional mechanisms may be involved, for example, 
FOXP2 (expressed at high levels in BP, but not control neurons) 
binds the DISCI promoter, 65 thereby potentially affecting cell 
migrations and cortical layering. There were no obvious differ- 
ences in the expression of transcripts for Hedgehog, Wnt or BMP 
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family members between groups, but the significant increase in 
transcripts encoding transforming growth factor-(3 family mem- 
bers in control iPSC may suggest a mechanism underlying the 
observed differences in dorsal/ventral neuronal phenotypes 
between groups. 

Immunohistochemistry demonstrated an increase in the num- 
ber of BP neurons expressing the MGE marker NKX2-1, whereas 
EMX2 was expressed by more control neurons. iPSC remained 
sensitive to patterning factors as demonstrated by treatment with 
the Hedgehog pathway activator purmorphamine or to LiCI. 
Interestingly, preliminary data also demonstrate that this early 
patterning is maintained with differentiation, with BP neurons 
expressing higher levels of somatostatin than controls. Consider- 
able research remains to be done to determine the timing and 
extent to which these alterations are maintained. 

The observation that DICER1 was significantly increased in BP 
neurons suggests a potential mechanism to explain the growing 
consensus that multiple genes may be misregulated in BP. As 
microRNAs alter the translation of as many as 70% of all human 
genes, 66 alterations in microRNA production would affect multiple 
gene products. MicroRNA levels are altered in postmortem brain 
tissue of BP patients, 67 ' 68 and are sensitive to lithium. 69 " 72 
MicroRNAs have critical roles in neural development, 73-75 as well 
as in the adult brain; 76 misregulation would therefore be expected 
to produce both short- and longer term affects on CNS 
development and function. 

Current work is in progress to expand the patient sample, to 
examine microRNA expression, to determine the extent to which 
patterning can be influenced by growth factors 42,43 or by the 
complex three-dimensional microenvironment of the differentiat- 
ing and mature brain, with the longer term goals of developing 
new animal models, and ultimately to identify new therapeutic 
approaches to this intractable disorder. 
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